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Abstract. Tiie decaying vacuum model (DV), treating dark energy as a varying vacuum, has been studied 
well recently. The vacuum energy decays linearly with the Hubble parameter in the late-times, pA{t) oc H{t), 
and produces the additional matter component. We constrain the parameters of the DV model using the 
recent data-sets from supernovae, gamma-ray bursts, baryon acoustic oscillations, CMB, the Hubble rate 
and x-rays in galaxy clusters. It is found that the best fit of matter density contrast Orn in the DV model 
is much lager than that in ylCDM model. We give the confidence contours in the Qm — h plane up to 3a 
confidence level. Besides, the normalized likelihoods of Om. and h are presented, respectively. 

PACS. 95.36.-t-x, 98.80.Cq, 98.80.Es 



1 Introduction 

Since the observations of the luminosity-redshift relation 
dL{z) of type la supernovae (SN la) indicated that the ex- 
pansion of the universe is accelerating [T, a cosmic com- 
ponent called dark energy was introduced to explain the 
acceleration within the framework of general relativity. 
Whereafter, more and more evidences, such as cosmic mi- 
crowave background (CMB) [5], large scale structure [5] 
especially baryon acoustic oscillations (BAO) [4], weak 
gravitational lensing [5] and x-ray clusters [Sll^ i indicated 
that the universe is spatially flat and dominated by dark 
energy. A number of dark energy models have been pro- 
posed [5], including scalar field |9 , vector field [TUIITT] . 
holographic dark energy [T2] , Chaplygin gas [13] and so on. 
Among them the cosmological constant model (ylCDM) 
[H] is the simplest one. However, as well known, the ylCDM 
suffers from the fine tuning problem: the observed vac- 
uum energy density of order ^ 10^'*^ GeV^ is about 10^^"'^ 
orders of magnitude smaller than the value expected by 
quantum field theory for a cut-off scale being the Plank 
scale, and is still about 10'*''^ orders smaller even for a cut- 
off scale being the QCD scale [S] . Apart from dark energy 
models, modified gravity [ISllTB] can also explain the ac- 
celeration. However, we just focus on dark energy model 
in this paper. 

As an extension to ylCDM, the the decaying vacuum 
(DV) dark energy model was proposed [TTlfTB] , based on 



expansion of the universe leading to an additional pro- 
duction of the matter component. In the late-time with a 
quasi-de Sitter background, the vacuum density is pro- 
portional to the Hubble rate, pA{t) oc H{t). However, 
the equation of state for the vacuum is a constant value 
w — PA{t)/ PA{t) — —1, the same as that in the ylCDM 
model. Moreover, as a interesting feature, the late-time 
dynamics of the DV model is similar to ACT)M [iTllTS] . 



The DV model has been tested by observational data 
of SN la [12], and a joint data from SN la, BAO, and 
CMB [201I2T] . Moreover, the quasar APM 08279-^5255 
at 2; = 3.91 was used to examine the DV model |22) . 
and it was found that the DV model can greatly allevi- 
ate the high redshift cosmic age problem existing in the 
ylCDM model. In order to distinguish the DV model from 
other dark energy models at the late-time Universe, the 
statefinder and Om diagnostics of the DV model were also 
presented in [22 . In order to constrain the DV model as 
completely as possible, besides the observational data of 
SN la of Union2 [53], BAO [15] and CMB [25] used in 
the recent work [21] , we added the data from Gamma-ray 
bursts (GRBs) at z > 1.4 [23, Hubble rate [27'.'2g and x- 
rays in galaxy clusters [7] . The confidence contours in the 
fim — h plane with the best fit values of the parameters in 
the DV model will be given, using various combinations of 
the data. Furthermore, one-dimension likelihoods of i7„ 



the incomplete quantum field theory in the curved 4-dimensiognd h, respectively, are also analyzed with the correspond- 
space-time. In this model, dark energy is described by the errors up to 3a 
varying vacuum, whose energy density decays with the 
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In the following, we use the unit SttG = c = 1. 
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2 The decaying vacuum model 



In a spatially flat Friedniann-Lemaitre- Robertson- Walker 
(FLRW) space-time, the Friedmann equation is 



(1) 



where the total energy density pT = Pm + Pr + Pde is com- 
posed with matter, radiation and dark energy. H = a/a 
is the Hubble rate of expansion, where an overdot means 
taking derivative with respect to t, the cosmic time. The 
dynamic evolution of the total energy density is deter- 
mined by the following conservation equation 



PT + 3H{pT +Pt) = 0, 



(2) 



where pt stands for the total pressure. For the present 
epoch, the radiation component can be ignored, and then 
the Universe only consists of matter (baryons and dark 
matter) and dark energy. So pr = Pm + PA and pt — Pa- 
In the DV model, the vacuum density varies with time. 
Pa — A{t)^ and moreover, the covarience of the Einstein 
equation demands the equation of state w — so pA — 
—A{t). Thus, Eq. (0) can be written as 



Pm + 3i/p„ 



(3) 



showing that the decaying vacuum density A{t) plays the 
role of a source for matter production. If A is a constant, 
Eq.Q reduces to the evolution equation for matter in 
ylCDM model. To proceed, one considers the late-time 
ansatz for DV model PTIITS] : 



A = (tH. 



(4) 



where cr is a positive constant. According to Eqs. ([TJ 
and ([4]), cr is determined as a — SilAHo, where f^A is 
the present value of the relative dark energy density and 
Hq = 100 h km s~"'^Mpc~^ is the Hubble constant. Com- 
bining Eq. ^ and Eq. ([4]), one can easily get the evolution 
equation of the Hubble parameter 



2H + -aH = 0, 



which can be rewritten as 



2H' + 3H -a = 0, 



(5) 



(6) 



with a prime meaning ' = d/d{hia). Eq. ^ can be easily 
resolved to get a solution of the Hubble parameter [T7] 



-3/2 



(7) 



Therefore, the matter energy density p„i = — aH and 
the vacuum energy density A = aH have the following 
explicit forms 



Pm{a) 



nla^ 



Aia) ^ - 



3/2>3 3i2^a3/2' 



(8) 
(9) 
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Fig. 1. The dynamic evolutions of the energy density of matter 
and vacuum in DV model and /ICDM model, respectively, for 
a fixed value of the present matter density contrast Orn = 0.3. 



where we have chosen a(to) = 1. The first terms in Eqs. 
^ and ^ give the usual scaling of the matter and the 
vacuum, respectively, while the second terms describe the 
matter production caused by the decaying vacuum. Eqs. 
^ and (O make an smooth transition between the matter 
and vacuum epochs, which are illustrated in Fig. [T] as an 
example with — 0.3. For comparison, we also plot 
the cases for ylCDM. One can see that the equal matter- 
vacuum dominated epoch in DV model is earlier {Az ~ 
0.4) than that in ylCDM model. From another point of 
view, the evolution of the Hubble parameter in Eq. ([7]) as 
a function of the redshift can be rewritten as 



H{z) =Ho[l- Qm + ^rail + zf'^] 



(10) 



Note that, this expression is rather different from that in 
the standard ylCDM model, due to the matter production. 
In particular, if yl = and Qra = 1, one has H{z) — 
Hq{1 + z)^/^ leading to pm = 3i?o(l + z)^, as expected 
for the Einstein-de Sitter scenario. On the other hand, if 
Ha = 1 and f^m = 0, we obtain H{z) = Hq and A — 3Hq, 
i.e., the Universe is described by the exact de Sitter space- 
time and the dark energy density will not vary. In the 
following, we will discuss the constraints on the DV model 
using various observational data. 



3 Observational constraints 

In the following, we use the observational data from SN 
la, BAO, CMB, GRB, Hubble parameter and x-ray in 
the galaxy clusters to constrain the two free parameters 
(]7m, Hq) of the DV model. As the likelihood function re- 
lates to function as £ oc exp (— x^/2), one should min- 
imize the function to obtain the best fit values of the 
free parameters. 
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3.1 SN la 

As did in many literatures, we compare the observed dis- 
tance modulus Hobsizi) of the 557 SN la assembled in 
the Union2 compilation [53] with the theoretical distance 
modulus 

fJ-thiz\ps) :^m-M = 51ogio + 25, (11) 

where m and M are the apparent and absolute magni- 
tudes, respectively, the complete set of parameters is Ps = 
{^m,Ho), and 



dL{z) = {l + z) 



dz' 



(12) 



is the luminosity distance in a spatially flat Universe. For 
the SN la data, the best fit to the set of parameters ps 
can be obtained by minimizing the function 



557 



XsNiPs)^}^ -2 , (13) 



where Ui stands for the la uncertainty associated with the 
ith data point. 



3.2 GRB 

Gamma-ray bursts (GRB) have been proposed to be a 
complementary probe to SN la Their high energy 
photons in the gamma-ray band are most intense explo- 
sions immune to dust extinction compared to supernovae. 
So far, there are many GRB observed in the range of 
0.1 < z < 8.1, and the maximum redshift of GRB is ex- 
pected to be 10 or even larger [30]. However, there is a 
circularity problem in the direct use of GRB 151 , mainly 
due to the lack of a set of low-redshift GRB at z < 0.1 
which is cosmology-independent. Not only some statisti- 
cal methods have been proposed to alleviate the circu- 
larity problem |32) . but also the cosmology-independent 
method needing a large amount of SN la at z < 1.4 to 
calibrate GRB was proposed in |33] . Due to both the sig- 
nificant improvements in GRB and SN la, the calibration 
of GRB was updated [231 , which can be used to constrain 
cosmological models without the circularity problem. The 
59 calibrated high-redshift (z > 1.4) GRB were obtained 
with the observed redshift z and modulus p, is shown in 
Table 2 of Ref. [H]. Hence, the Eq. ^ is also valid for 
GRB, 



XGRsiPs) 



59 

E 

i=l 



[^J■thizi) - Pobs{zi)Y^ 



(14) 



3.3 BAO 

In the distribution of SDSS luminous red galaxies [MIH] 
the distance parameter A for the measurement of the BAO 



peak has been usually used as a kind of examination in 
literatures. The definition of A is 



A = Dx 



where Dv is the dilation scale, defined as 



Dv{z) 



H{z) 



H{z') 



1/3 



(15) 



(16) 



However, as pointed out in 20 , the distance parameter 
A is not appropriate to test the DV model, since it is 
obtained from the data in the context of a ylCDM model 
and can be considered as a good approximation only for 
some class of dark energy model [3F . We follow Carneiro 
et al. ^SOj, using Dy instead which is weakly sensitive to 
the cosmological evolution before z = 0.35. The dilation 
scale Dv has been observed by the SDSS at z = 0.35 
m, as weh as by the 2dFGRS at z = 0.2 [IS]. Moreover, 
the measurements of SDSS and 2dFGRS yield the ratio 
Dv{Q.ih)lDv{Q.2) = 1.736 ±0.065 Then the best fit 
values for the model can be obtained by minimizing |36) 



x]3AoiPs) 



([J„(0.35)/J„(0.2)]tfe ~ [A,(0.35)/J.„(0.2)]ob,)^ 

2 

'^_D^C0.35)/D^,{0.2) 



(17) 

where (t|, zb)/D (o 2) ~ 0.065. However, as presented in 
Ref. [U, it does not impose any constraint on h. 



3.4 CMB 



The CMB shift parameter defined as [571155] 



dz 



(18) 



with Zrec the redshift of recombination, often serves as 
a contrast in a large class of dark energy models to be 
compared with that in the ylCDM model. However, this 
is only valid if the acoustic horizon at the time of last 
scattering is the same [33] ■ This is not true in DV model, as 
pointed in [20] , because for the same values of Hq and fim 
we have different expressions for cosmological parameters 
at high redshifts, due to the process of matter production. 
In this paper, to employ less fitting expressions, we use the 
acoustic scale I a instead of the position for the first peak li 
in the spectrum of CMB anisotropics [^[^ . The acoustic 
scale at decoupling epoch is defined as 



_ fz* dz 
J _ ^ Jo Hjz) 

Jz, H(z) 

where is the redshift of decoupling epoch, and 



(19) 



4 Q^{1 + z) 



-1/2 



(20) 
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is the sound speed of the photon-baryon fluid. Here J7b and 
Q~f are the present relative energy densities of baryons and 
photons. Below, we follow the observational results in [26 : 
= 302.09, = 1091, Db = 0.0455 and = 2.469 x 
10~^/i"2 for TcMB = 2.725K. Moreover, at the very early 
stage for the DV model, the appropriate generalization of 
Eq. (ITOl) including radiation is given by [20] 



H{z) 



{[l-l2„ + l?™(l + z)3/2]2 + x2,(l + z)4}i/2^ (21) 



where f2r is the present relative energy density of radia- 
tion. Here, f^r = fi'y{l + 0.227 x Nef f) with the relative 
energy density of photons and -/Ve// =3.04 the effective 
number of the standard neutrino species |40| . Therefore, 
the corresponding function is given by 



XcMB^Ps) 



with cr;^ = 0.76 



{I 



A.th 



lA,obs)'^ 



(22) 



3.5 Hubble rate 

Recently, some high precision measurements constrained 
H{z) at z = from the observation of 240 Cepheid vari- 
ables of rather similar periods and metallicities [57| . The 
Hubble parameter as a function of redshift can be written 

as 

H{z) = --^^. (23) 
1 + z at 

Then H(z) is obtained once dz/dt is known. Simon et 
al. [H] and Stern et al. [42] obtained H{z) in the range 
of < z < 1.8, using the differential ages of passively- 
evolving galaxies and archival data. Besides, H{z) at z = 
0.24, 0.34 and 0.43 is obtained [28] by using the BAO 
peak position as a standard ruler in the radial direction. 
We employ the twelve data in |?flH^ and the three data 
in [28]. The best fit values of the model parameters from 
observational Hubble data are determined by minimizing 



15 

XHubiPs) = 

1=1 



[Hth{z^) - Hobs{zi)? 



cJ^zi) 



(24) 



3.6 X-rays in galaxy clusters 



The baryons in clusters of galaxies are in the form of hot 
x-ray emitting gas clouds. Thus, the fraction of baryons 
in clusters, fgas, is defined as the ratio of the x-ray gas 
mass to the total mass of a gas mass fraction, which is a 
constant and independent of redshift [33]. Here, we use 42 
Chandra measurements of relaxed galaxy clusters in the 
redshift range 0.05 < z < 0.1 [J. To fit the data one can 
employ the empirical formula |36| 



fgas{z) 



KA-fb{z) 

1 + s(z) 



V dA 



3/2 



(25) 



where the angular correction factor A is approximately 
given by 

H{z)dA{z) 



[H{z)dA{2 



]ACDM 



(26) 



Here dA is the angular diameter distance, being related 
with the luminosity distance by 



dAiz) 



[l + zf 



(27) 



and the index 77 = 0.214 ± 0.022 f7]. In equation (^5)) . 
K describes the combined effects of the residual uncer- 
tainties such as the instrumental calibration and certain 
x-ray modelling issues. The parameter 7 denotes permis- 
sible departures from the assumption of hydrostatic equi- 
librium due to non-thermal pressure support. The bias 
b{z) = &o(l + abz) accounts for uncertainties in the clus- 
ter depletion factor, and s(z)so(l -|- agz) answers for un- 
certainties of the baryonic mass fraction in stars. Most of 
these parameters have large uncertainties which can be 
seen in Table 4 in |:7j. For simplicity, we fix these param- 
eters in the following: K = ^ = 1, sq = QA&h^/"^, as = 0, 
&o = 0.835 ^i"', and ab — 0. The function from x-rays 
in galaxy clusters is then given by 



siPs) 



42 

E 



( \_f gas (^-i)] \_f gas {,^i)\ohs ) 



(28) 



4 Results and discussions 

Based on the discussion above, the total is combined 
as 



A^total 



XSN 



2 

Xgrb 



2 I 2 I 2 

XbAO + XCMB + XHub 



Xx 



-rays ' 

(29) 



As the likelihood function is determined as £ oc exp(— x2/2), 
the left panel of FigH shows the 68.3% (Icr), 95.4% (2a) 
and 99.7% (3fT) confidence contours of the parameters 
(f^m, h) obtained by various combinations of observational 
data-sets for the DV model. In order to find how much 
a set of data contributes to the constraints, we added 
the sets of data one by one. We can see that the data 
from CMB and x-rays give very sensitive contributions to 
the constraints. Here we should emphasize two aspects. 
Firstly, as discussed in |20| , we have used parameters that 
are strictly correct for the ylCDM case, which would lead 
to bias in our results in spite of a good approximation. 
So, more complete analyses of CMB need to be done. Sec- 
ondly, as can be seen in Eq. (^5]) . the value of Xx-rays 
depends on the ratio of the luminosity distances predicted 
by ylCDM and ESF, which will deviate from 1 more and 
more at higher redshift. This is why the result of the best 
fits is so different with and without the x-rays data. Fur- 
thermore, there are many uncertain parameters included 
in Eq. (|25t . Different choices of the parameters lead to 
different contributions from x-ray data to the results. On 
the other hand, the additions of the data from GRB and 
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Fig. 2. Left panel: The confidence contours up to 3a in ilm — h plane for various combinations of observational data. The 
constraints from total data are shown as shaded contours. Right panel: the constraints for the DV model and /ICDM model 
from the total data and the combined data without x-rays, respectively. The solid points stand for the location of the best fit 
values. 



Hubble parameter give negligible influences on the con- 
straints. One can see from Fig.[5]tfiat, the red dash-dotted 
contours and tlie black dotted contours overlap with each 
other, and the green solid contours almost overlap with 
the previous two sets of contours. Concretely, the best fit 
values of parameters f^m and h with corresponding Xmin 
and Xmin/^'^f ('^o/ = degrees of freedom) obtained from 
different combinations of data-sets are presented in Ta- 
ble [TJ respectively. In Table [U we can see that the best 
fit values of f2rn and h are very sensitive to the data from 
CMB and x-rays, but are very insensitive to the data from 
GRB and Hubble parameter. For comparison, we also give 
the two cases in the ylCDM model. A comparison of confi- 
dence contours between DV and ylCDM can be seen in the 
right panel of Fig.[2J The lined contours are obtained from 
joint data of SN+BAO+CMB-f GRB-|-Hub, and shaded 
contours are obtained from the total observational data. 
The results of the confidence contours obtained without 
x-rays data, as can be seen in the right panel of Fig. [5] are 
in good agreement with those appeared in [201I2T] . 

Fig. [3] presents the 1-dimension normalized likelihoods 
of the two parameters i?^ and h for DV model based on 
the total data (blue lines) and the combined data without 
x-rays (red lines) , respectively. Concretely, the best fit val- 
ues are = 0.3960t°:°^34(i^)+o.oo70(2a)t« ojo4(3^) 

h = 0.6978±0.0025±0.0050±0.0074withxLai = 966.343. 
If we abandon the data from x-ray, the best fit values 

arp O - n 4970+0 0044+0. 0089+0. 0132 u - Q CQco + 

are j/^ — u.4z^u_q 0043-0. ooss-o. 0125 " — u.osos m 
0.0025 ± 0.0050t[J:;]J]^^ with Xtotai = 596.345. One can see 
the big different results of the both likelihoods of Sim and 
h. As discussed above, the differences are due to the devi- 
ations of the luminosity distance given by ESF from that 
given by ylCDM. 

Just as the results from the examination using en- 
ergy density perturbations [JS], the best fits of J7m for 
DV model are very large compared to those for ylCDM 



model. The data from CMB give a magnificent response 
for the large value of However, more complete anal- 
yses of CMB need to be done. Moreover, there are many 
parameters with large uncertainties in Eq. (I25p . when we 
evaluate Xx-rays in Eq. (pS)) . This defect can just explain 
the large value of x^^^/dof = 1.44, which can be seen in 
Table ffl 



5 Conclusions 

We constrained the parameters of the DV model using 
the most recent joint data from 557 SN la in Union2 
compilation, 59 GRBs in the high redshifts {z > 1.4), 
the ratio of the dilation scale at z = 0.35 and z = 0.2 
of BAO, the acoustic scale at decoupling epoch of CMB 
measured by WMAP 7, 15 data of Hubble rate, and 42 x- 
ray gas mass ratio in galaxy clusters. Due to the joint 
data, the parameters of the best fit values are: flm = 



0.3960 



+0.0034+0.0070+0.0104 



and/i = 0.6978±0.0025±0.0050± 



-0.0033-0.0067-0.0099 

0.0074. With the same joint data, we also calculated the 
best fits for ^CDM model: = 0.293lt° ™i^ll° ™^?t° °J°^ 
and h = 0.6941 ± 0.0024 ± 0.0047 ± 0.0070. However, the 
data from x-rays depend on many fitting parameters with 
a large uncertainty. Then, if we neglect the poor data of 

X ravs the results are- O - 4070+0 0044+0. 0089+0. 0132 

X-rays, tne results are. Hm — u.^z (^u_o 0043-0. 0085-0. 0125 
and h = 0.6858 ± 0.0025 ±0.0050t;^:|^^^^ for the DV model; 

n„, = 0.2785l^ ™^^l!; °ii^^tni^t and h = 0.6983±0.0027± 
0.0053 ± 0.0079 for the ylCDM model. The best fit value 
of J7m in the DV model is larger than that in the ylCDM 
model, since there is additional production of the mat- 
ter component due to the decaying of the vacuum energy. 
However, we cannot rule out the DV model from the cur- 
rent observational data. We should wait for more abun- 
dant, accurate, and model-independent data to constrain 
this model. On the other hand, the DV model alleviates 
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Table 1. The best fit values of the DV model using different combinations of data-sets. 



Data 




h 






SN 


0.3388 


0.6953 


544.587 


0.98 


SN+BAO 


0.3369 


0.6955 


545.758 


0.98 


SN+BAO+CMB 


0.4272 


0.6857 


561.724 


1.01 


SN+BAO+CMB+GRB 


0.4272 


0.6857 


584.866 


0.95 


SN+BAO+CMB+GRB+Hub 


0.4271 


0.6859 


596.345 


0.95 


total 


0.3960 


0.6978 


966.343 


1.44 


SN+BAO+CMB+GRB+Hub (ylCDM) 


0.2785 


0.6983 


577.915 


0.92 


total (ylCDM) 


0.2931 


0.6941 


632.412 


0.94 




Fig. 3. The normalized likelihoods of the parameters and h, respectively, for the DV model. 



the fine tuning problem and high redshift cosmic age prob- 
lem which exist in ylCDM model. 
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